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On the Determination of a certain Glass of Inequalities in the Moon's 
Motion. By Ernest W. Brown, B.A. 

1. The two principal methods of treating the Lunar Theory— 
viz. (1) the general method, in which we obtain a general ap¬ 
proximate solution, applying to any single Moon whose present 
orbit lies within given limits, in terms of certain constants, the 
coordinates of position and the time; and (2) the special method, 
in which we use the numerical values of those constants at the 
outset for any given Moon, and obtain a solution involving the 
coordinates of position and the time only—present each of them 
certain advantages. In the first case, we have the theory, worked 
out to a certain degree of accuracy, immediately applicable to 
any single Moon in our solar system, and therefore arranged in 
such a way that any small change which improved data may 
involve in the values of the constants can be made easily with¬ 
out requiring us to go over the whole of the work again. On 
the other hand, the number of terms which have been found 
necessary to secure a degree of accuracy commensurate with that 
of observation is very large, and it becomes a task of great labour 
to obtain them with any degree of certainty. In the second 
case, when we use numerical values from the start, the portions 
arising from the expansion of the functions in ascending powers 
of the constants involved, which portions in the first case would 
have to be neglected, are no longer left out of account. A great 
increase of accuracy naturally results. But this method is un¬ 
satisfactory in some respects. Having started with certain values 
of the constants for any particular Moon, we are bound to keep to 
those values ; and also if any numerical mistake be made, it is not 
very easily traceable. The calculation is applicable only to the 
particular Moon selected. Any method, therefore, of finding the 
path of the Moon which will give at the same time numerical 
and algebraical results with equal accuracy would be of value. 

Mr. G. W. Hill, in the American Journal of Mathematics , 
vol. i.,* has shown how to separate out those inequalities de¬ 
pendent only on the ratio of the mean motions of the Sun and 
Moon, and to obtain their coefficients b6th numerically and as an 
algebraic function in ascending powers of this ratio. He shows 
also that his method is not only susceptible of a very high degree 
of accuracy with comparatively little trouble, but that any higher 
degree of accuracy can be obtained without the necessity of going 
over a large part of the ground again. His method for this 
class of inequalities is as follows:— 

2. Take as origin the centre of the Earth considered as a 
sphere, and refer the motion to rectangular axes revolving round 
the origin with the mean angular velocity n! of the Sun round 
the Earth. If we neglect all the differences from purely elliptic 

* 'Researches in the Lunar Theory , pp. 5, 129, 245 
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motion, but those dependent on the mean motions of the Sun and 
Moon, the Sun will describe a circle about the Earth with con¬ 
stant angular velocity n ! , and we may therefore assume the 
moving axis of x to pass through the Sun. Under these circum¬ 
stances, the equations of motions of the Moon become 


dfa 
dt 2 


— 2 n f ^ + 
dt 


(£-*■>- 


d?y 
dt 2 


+ 27l f 


dx u. 


= 0 


I. 


These equations will give us all inequalities from purely circular 
motion, which involve only the eccentricity of the Moon and the 
mean motions of the Sun and Moon. By taking certain “ par¬ 
ticular integrals” of these equations, the inequalities dependent 
only on the mean motions can be separated out and determined. 
This particular solution is given by 

x = '%a 2i cos ( 2 z + 1) (n — n ') (t — 1 0 ) 
y = sin (21 + 1) [n — n’) (t — 1 0 ), 

when i takes integral values positive and negative, and a 2i is a 
function only of the ratio of the mean motions of the order m 2i 
at least. The quantities a 2i are then determined by substituting 
this particular solution in the equations of motion. Certain 
transformations are first made in order to reduce the equations (I) 
to a convenient form. 

3. In most previous methods, and especially in Delaunay’s, 
the coefficients of this class of inequalities are expressed directly 
in ascending powers of n'/n=m; but in Mr. Hill’s method of 
preparing the equations for solution the function of m which 
appears most naturally for use is n' / (n—n'f Let 

n' rn 

m, = -,=-; 

n — n i—m 

we then have 

m, 

m — - 1 . 

I + Wj x 

At first sight m l would seem to be a less suitable quantity than 
m, since is greater than m, and consequently a greater number 
of terms would seem to be required to obtain a given degree of 
accuracy. But it is shown that the numerical multipliers of the 
higher powers of m x are much smaller than those of the corre¬ 
sponding series in m ; in other words, the series for the co¬ 
efficients expressed in ascending powers of mj are so much more 
convergent than the same series expressed in powers of m, as to 
quite overbalance the disadvantage of m 1 being slightly greater 
than m. Further, he shows that by using a certain function of 
m x obtained theoretically , for expansion, the series can be ren- 
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dered still more convergent. This function of m l is /x, where /a 
is given by 

_ m l __ m 

i->' 


Hence 




i+l^’ 


The series for a 2 i are given by Mr. Hill, developed in terms of 
both m x and /x. 

4. These transformations can be utilised also in another way. 
The inequalities found as indicated in § (2) are referred to rect¬ 
angular coordinates. It has been usual to refer the Moon’s 
motion to polar.coordinates. We have for this purpose 

r cos v = Za^ cos 2 i (n—n') {t — tt 0 ) 
r sin v = Za^ sin 2 i (n—n')(t—\ f 0 ), 


where r is the radius vector and v the excess of the true over the 
mean longitude. From these r and v can be easily found. 

The process of finding these coefficients as algebraic expan¬ 
sions in powers is somewhat laborious when, as in the special 
case of the Variation Inequality, it is necessary to carry the 
expansions to a high power. But there is no reason why, 
supposing Delaunay’s expressions in longitude and parallax for 
these coefficients to be accurate as far as they go, we should not 
put in them 



and expand in powers of /x. This transformation gives us smaller 
numerical multipliers of the higher powers, and consequently 
greater accuracy when we substitute for /x its numerical value. 
It is to be noted that if the given series is calculated as far as 
m 2 , the transformed series cannot be carried further than /x 2 . 

5. As an example, take that part of Delaunay’s expression* 
for the coefficient of the Variation which depends on m. It is 


11. 

S 




+ ^m 3 
12 


+ ?23 ro 4 
72 


+ 28 55 m . 

IO8 


8304449 6 
165888 


I586"*8883 +424"*4474 + 8o"-o9o6 +12"7689 + i"-8o87 


1 102859909 ^ 
1244160 

+ o"-2234 


+ 759 6606 727^, 
74649600 

+ 0"-0206 


_8o5I4i8i6i 
1119744000 

— o n, oooi 




(a) 


where the coefficients expressed in seconds are written under¬ 
neath. Putting in this series 

m = nl(i+ &), 

* Delaunay, Memoires de VAcademie des Sciences , tome xxix. p. 815. 
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and expanding in powers of /z, it becomes 


ii 


+ S -A 


+ -^-/z 4 
72 


~^ 5 

36 


i957"'9686 + i47"-8944 + o"-6827 — o"-2496 


82111 6 

“ispsT 

-o"-0336 


_350399^7 
138240 


— o f, -oi43 


_233559113^8 
74649600 

- o n -ooi5 


_10961275281^ (0) 

1119744000 

— o f '*ooo4 


The values of tlie coefficients obtained from these expressions 
are— 

From (a) 

From (£) 

From Mr. Hill’s values 


n 

.. 2106-2478^ 

.. 2106-2463 1 .(7) 

.. 2106-2463] 


The last value has been calculated from the numerical values of 
the coefficients a 2i given in his paper already referred to. The 
equality of the second and third values—the second obtained 
from Delaunay’s series when expressed in terms of jjl, and the 
third from Mr. Hill’s numerical values—is striking. It not 
only shows us the improvement in accuracy afforded by such 
transformations, but it gives us also some security for the general 
correctness of this particular part of Delaunay’s expression for 
the coefficient of the variation. 

6 . There appears to be a kind of oscillation in the class of 
series of which (a) or ( J 3 ) is an example. For instance, in the 
series (a) the numerical multipliers are positive, and increase 
gradually until a change of sign occurs between m 8 and m 9 , 
when there is a sudden drop. The numerical values (y) seem 
to show that there must be a large increase in the multipliers of 
the powers after m 9 , and that they must have negative signs ; so 
that probably the multipliers of m 10 and m 11 are much larger 
than even those of m 7 or m 8 . The series (/ 3 ) has a sudden drop 
between /z 3 and /z 4 , and then a change of sign occurs, after which 
the multipliers gradually increase again. I have calculated 
similar expressions in the coefficients of other inequalities, and 
generally it appears that when we, by any such substitution as 
that made above, apparently improve the convergency of the 
series, the change of sign is brought nearer to the beginning of 
the series. An ideal to aim at would seem to be, firstly, that 
the numerical multipliers be made as small as possible; secondly, 
that there be no sudden increase in them in the later part of 
the series; and thirdly, that these two conditions should involve 
that there be no long run of powers with the same sign attached 
to them. This appears very forcibly in the expressions for the 
Parallactic Inequality given in § (i i) below. It should be stated 
that as the series are calculated up to some definite power only* 
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and as the law of progression of the series is not able to be 
expressed by an algebraical formula, owing to the complicated 
forms from which they arise, a substitution like that made above 
must not be arbitrary , but must be indicated by theory. An 
arbitrary substitution, which would seem to make the parts of 
the series already calculated converge rapidly, may make the 
unknown part very slowly convergent, and thus introduce un¬ 
known errors into the numerical values of the coefficients. 

7. It has been mentioned above that the equations (I.) deter¬ 

mine all the inequalities dependent on the eccentricity and the 
mean motion only. As in these equations when transformed for 
solution m x and not m is the quantity which appears, all those 
series which depend on the eccentricity and mean motions only 
will probably be rendered more convergent when expressed in 
terms of m x . Consequently, in Delaunay’s expressions involving 
e and m only, we are entitled to put +%), and expand 

in powers of m 1 up to the same power of m Y as that to which m 
was carried; the result of this should be more convergent series 
and greater accuracy in the coefficients when expressed in 
seconds. It will be found on examination that nearly all these 
series have a number of successive powers of m of the same sign 
with gradually increasing multipliers; and the transformation 
improves the series in the manner mentioned in § (6). 

8. When we examine Delaunay’s expression for the part of 
the coefficient of the Parallactic Inequality dependent on the 
ratios of the mean distances and mean motions, we notice that 
the numerical multipliers of the higher powers of m become very 
large. The multiplier of any one power is between four and 
seven times that of the previous power; and since m is in the 
case of our Moon about *074, the series expressed in seconds 
converges very slowly. The expression is given below in § (11). 
In fact, the last term calculated, which is of the ninth order, 
corresponds to a coefficient of o"'38; so that some method is 
necessary for this inequality which will give us either more 
terms easily, or a new form of expression with greater con- 
vergency. I shall give elsewhere the analysis by which the 
expressions and values below have been obtained, and expose 
only the method used to work up the coefficients into a more 
convergent form. No other inequalities but those given are 
affected thereby. 

In obtaining these inequalities dependent on m and the ratio 
of the mean distances, some terms of the disturbing function 
have to be added to the left-hand sides of the equations (I.). The 
constant m 1 is still the function of m which appears naturally 
for use in expansion. The general method of procedure has 
been the same as that Mr. Hill used in getting the inequalities 
dependent on the mean motions only. 

9. Let a denote the mean semi-axis major of the Moon’s orbit, 
a' that of the Sun’s orbit, v the difference between the true and 
mean longitudes of the Moon, 2D the argument of the variation 
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in Delaunay’s theory; the parts of the Lunar Inequalities 
dependent on the mean motions and a /V only, appear in the 
following form:— 

r cos v — a 0 + (a l + a-i) cos D + (a 2 + ti_ 2 ) cos 2 D + ... 
rsint;= «._j) sin D + (a 2 — a_ 2 ) sin 2 D + ... 


The quantities a± v are functions of the kind, 

where jp is a positive integer. 

In the functions / for a 2p , m x p is the lowest power of m l 
which occurs; in F the lowest power is also m x 2p , except 
when p=o> when it is ; when p=o, f is of the form 
1 +am 1 2 + / 3 m 1 3 + . . , a, . . being numbers. Mr. Hill, in his 
paper referred to, has given the functions /, and his values for 
them have been used in obtaining the functions F, <£, <£>. 

When we transform these expressions to polar coordinates, 
we have, neglecting cubes and higher powers of &/&', 

v = A ~ sin D+ |b + B‘ j s i n 2D + C^- r s ^ n 3L + . • . 

where A, B, B ; , C . . . are functions of m 1( D is the argument 
of the Parallactic Inequality, and the first term represents the 
principal part of it, i.e. the part dependent on m x and aja' only. 
The full expression of this inequality contains terms dependent 
on the solar and lunar eccentricity and the lunar inclination, but 
we are not here concerned with such terms, and they do not 
affect the expressions which follow. 

10. The principal part of the coefficient A aja' depends 
chiefly on a x /a and a^ja. In obtaining a Y and a_ x as far as the 
order m x 5 . aja', it is possible so to arrange the work that the 
equations determining them are linear and of the form 


Cb 1 CL _i « 

1 . «,+ —1 • 0 !=*! 
a 0 


a 

a' 


a , a-, 

- 1 - . a 2 + _ 1 
a Q a 0 


— K 2 


a 

a’ 


where a T a 2 fix K 1 K 2 are functions of m 1 expanded in ascend¬ 
ing powers of m t , and the signs in the equations are so arranged 
that these functions become positive numerical quantities in the 
special case of our Moon. 

Solving these equations for a x and a_ x we have 

a x _ _/c 1 j8 2 + /cjS 1 a a_ 1 _ K x a 2 +K 2 a x 

Cb 0 ^i a 2 — 0 lAs O' a* ^ 
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The expressions on the left-hand sides of these equations, 
when for a Y /? 1? etc. are substituted their expansions in powers 
of are cumbrous. At the same time, when a 1 fix, etc. are 
given their numerical values, we get a very close approximation 
to the results obtained by the more accurate method of using 
numerical values from the start. Each of these quantities, when 
expanded in powers of m l3 converges quickly and no very large 
numerical multipliers seem to appear. But the multipliers in all 
of them gradually increase, and consequently, when we express 
the products /q /? 2 , *2 etc., i n ascending powers of m l9 larger 
multipliers occur. It seems best, however, for simplicity to 
expand the numerators and denominators of the expressions for 
aifa Q and a-i/a Q , and leave the results in the form of fractions 
which have the same denominator. The loss in accuracy is very 
small. We have then 

a x _ _ 7 (^i) a a ~\_ + ^( w 0 a 

a 0 \ r a' a 0 r * d 


where y, S, r are functions of m 1 expanded in ascending powers. 

11. The terms in longitude with arguments D, 3D, 5D, 
calculated in this way, become 


- *5, 


-m, +z.m 1 z — _„ 

8 1 4 1 I2B 


*95^ 3_ 41585™ 4 _ 2096751 


1024 


■mr —. 


49152 


- m , 


1 —4§ 9963 4 
o 6 9216 


■ - 

1024 12288 


■i24 /, -i493 -i2"*0447 + 6"'597o + i"'42i8 

-0"'0008 — o"'CKX)2 


+ I5 m 2 + 55 m s_ _4i m 4 + £3^9 w 5_ 
32 1 12? 1536 9216 


x ^Lsin D 
a' 


+ o”i207 


255 s + 7543 4 _ I345 2 5 5 
128 1536 8192 

i-4^i . . . 


x -L sin 3D 
a’ 


+ i"-6i 72 +0-1199 -o"ooo6 +o"-ooo5 

— o”-8624 -o"i7i8 + o"*047o 


■ 1W + 27i§4is^. x _1_ 

+ 128 2048 8 1 92 . . . 

+ o ,f -oo93 


ci • Tv 

x — sm 5D. 
x a r D 
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For the sake of comparison I give Delaunay’s expressions* for 
the same terms; they are expanded in powers of m. 

- + + 6887 m3 + WW*,. + 462833 3^ + 63106813^ 

[_8 8 128 512 3072 8192 

-74"-0235 - 34"'3297 -ii"* 8852 -4 r '*4276 -i"‘862i -o"*7i22 

+ 10835537159 m 7lg_ gin T) 
196608 -J «' 

-o"-38ii 

+ fi^m 2 -|m 3 -?^m 4 -^? 49 m 5 l-sin 3D 
L32 8 16 6144 Jo' 

+ i"‘3843 —o ,, *i38i — o"* 2675 — o"*i340 

+ s i n 

L128 1024 J a' D 

+ o"-oo97 + o''*0047 


The gain in accuracy by the former mode of expression is 
evident; the terms there are carried only as far as Wj 5 , while 
Delaunay’s go as far as m 7 . 

The equations for obtaining oq involve largely a 2 a ^ 2 , 
which go to form the coefficient of the. Variation. It has been 
mentioned that the expansions of these quantities are much more 
convergent when expressed in terms of //,, § (3). Hence it is 
probable that the resulting expressions for the coefficients of the 
Parallactic Inequality would be also rendered more convergent 
when expressed in terms of this quantity. This inequality thus 
transformed becomes, 


is + _ 2I| M ,_ 175705 4 + 3261161 , 

8 ^_128_9216_442368 ^ 

I _^_i33^ + 9I.^_7oi89 i , 

3 72 216 82944 


+ ^ 35 _ fJL ± + 
IO24 


? 8 5 


12288 


CL • t\ 

_sm D. 
a' 


— I3i' /, i3ii _^"- 8 io 9 + 8 "-i 438 +o "*7648 — o "-0246 — o"*ooo 9 — o"*oooi 


[In making the transformation, both numerator and denomi¬ 
nator of the expression given at the beginning of this section 
must be multiplied by m 1} and, in the result, for m l must be put 
/ut/ (1 + 37*). It is i n t^ a t form that the expression occurs.] The 
series in both numerator and denominator seem slightly more 
convergent, and the numerical result is a little closer to the final 
value. 

12. To exhibit the values of the coefficients obtained by using 
these several expressions, a table is given. In column (i.) is the 
argument of the term ; in (ii.) the coefficients found by using 
the numerical value of m x from the start; (iii.) those found by 

* Delaunay, M 6 m . cit . pp. 847, 856, 860. 
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putting the numerical values of a Y /? 1? etc. in the expressions of 
§ (10) ; (iv.) gives the values from the expressions in terms of 
m l in§ (ii.); (v.) those in terms of in the same section ; in (vi.) 
are Delaunay’s values. No estimate of neglected terms is taken 
into account. 


i. 

ii. 

iii. 

iv. 

V. 

vi. 

D 

— 128-069 

- 128-073 

— 128-056 

-128-059 

/ 

— 127621 

3D 

+ 0*750 

... 

f 0-749 

... 

$ 

b 

+ 

5D 

+ 0-008 

... 

+ 0*009 

... 

+ 0-014 


The column (ii.) gives the final values to be used; (iii.), (iv.), or 
(v.) may serve also as checks on the numerical work. It may 
be noticed that in obtaining the numerical values from the outset 
it is only necessary to substitute for m l its value. The ratio a/a' 
occurs only in its first power; it may therefore be kept as a 
factor right through, and its numerical value substituted as the 
final step. 

13. The expansions I have given for the parallactic terms, 
when transformed, agree with those of Delaunay as far as the 
terms of the sixth order. The numerical multipliers of 

differ by small quantities; these differences would not, however, 
cause differences in the coefficients expressed in seconds of more 
than o // *oo4. As the coefficient of the term with argument 5D 
is very easily obtainable by the methods indicated, I am inclined 
to attribute some small numerical errors to Delaunay. Newcomb 
suggests * that Delaunay’s multiplier of m 7 u/a' may be wrong. 
I have not yet been able to test Delaunay’s multipliers of m 6 . a/a! 
and m 7 .a/a f owing to the large amount of labour they entail. 
But that the numerical values given in column (ii.) of the table 
do not differ by more than one or two thousandths of a second 
from the true values can be seen during the progress of the work 
necessary for finding those values. 

14. The corresponding terms in the expression for the Moon’s 
parallax, which I have obtained, are :— 


1 

a 


i 1 ? 3 o 1869 

_s'/m _L_ m 777, * _ 

256 


1r* + s^ 




I — 4 ?n — . . . 


CL tv 

-j cos D. 


— i"-029i -o''’0333 +o"*o524 = — i"-oioo. 

+ - Rr* m l 2+ — . -—-COS 3D. 

a L64 256 256 1 — . . . J a 

+ o"Q223 +o"-ooi2 —o"-oi39 =+o"-oo96. 


* Newcomb, Transformation of Hansen's Lunar Theory , Astronomical 
Papers for Use of American Ephemeris, vol. i., part ii. 
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Delaunay’s expressions * for the same being 

i TK 8i a , 5817 „“1 a -p. 

— ~^-m + — m'? + ~ — L m 3 — cos D. 

a\_i 6 l6 256 J a' 

— o”*6i36 -o"'2479 -o"*o 832 = —o''*9447. 

+- f— ni 1 — iiirn 3 ! - cos 3D. 
a L64 128 _U' ^ 

+ o"-oi9i —o"-oo33 = +o”*oi58* 


LII. 2, 


The coefficients obtained by nsing numerical values from the 
start are — i'^ooio and +o // *oo8o, to which the former pair 
approximate closely. 

15. A comparison of these new values with those given by 
Hansen cannot be complete owing to the fact that Hansen’s 
numerical values contain the terms due to eccentricity and incli¬ 
nation. But a rough estimate of the uncalculated parts of these 
series may be made from Delaunay’s expressions, since they are 
carried to a higher degree of approximation for numerical pur¬ 
poses than the terms we have been considering. The whole co¬ 
efficient of the Parallactic Inequality in the longitude then 
becomes 

i27"*3o 


Newcomb’s method of transformation f shows that this coeffi¬ 
cient when expressed in Hansen’s notation is 

124-54 

Hansen’s coefficient is 124*43 
Diff. o*n 


The difference is probably partly due to a wrong estimate of the 
omitted terms in Delaunay’s coefficient, 

Haverford Coll., Pa., U. 8 .A. 


Secular Perturbations of the Earth's Orbit by Mars. 

By R. T. A. Innes. 

In Professor Hill’s li Theory of Jupiter and Saturn ” he men¬ 
tions that the terms arising from the fifth powers of the eccen¬ 
tricities and inclination amount to one per cent, of those of the 
first order in the case of Mars and the Earth, and that a com¬ 
putation by Gauss’s method would be desirable. I have therefore 

* Delaunay, Mem. cit. pp. 922, 923. 
f Newcomb, Mem. cit. 
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